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Achieving controllable coupling of dopants in silicon is crucial for operating donor-based qubit 
devices, but it is difficult because of the small size of donor-bound electron wavefunctions. Here 
we report the characterization of a quantum dot coupled to a localized electronic state, and we 
present evidence of controllable coupling between the quantum dot and the localized state. A set of 
measurements of transport through this device enable the determination of the most likely location 
of the localized state, consistent with an electronically active impurity in the quantum well near the 
edge of the quantum dot. The experiments we report are consistent with a gate-volt age controllable 
tunnel coupling, which is an important building block for hybrid donor and gate-defined quantum 
dot devices. 


Donors in silicon are a natural choice for qubits 1 , be¬ 
cause electron spins bound to donors have very long 
coherence times 2-4 . Phosphorus donors also have nu¬ 
clear spins with extremely long coherence times 5 7 . Al¬ 
though donor-based quantum devices can be fabricated 
with near-atomically precise placement of donors 8,9 , even 
when well-placed, donors are very small, making it diffi¬ 
cult to control and change the tunnel couplings between 
them with gate voltages. In contrast, tunnel couplings 
are easily tunable in gate-defined quantum dots, and 
high-quality quantum dots hosting at least four differ¬ 
ent types of spin qubits have been demonstrated semi¬ 
conductor materials 10 22 . An important feature of gate- 
defined quantum dots is that the electrons they contain 
can be displaced laterally simply by changing the volt¬ 
ages of the gates on the surface 23 . Because of the ease 
of spatial control of the wavefunction in quantum dots, 
coupling quantum dots to localized states in semiconduc¬ 
tors, if possible, would enable a path of a wide range of 
hybrid donor/quantum dot technologies. 

In this letter, we report the observation of a control¬ 
lable tunnel coupling between a localized electronic state 
and a gate-defined quantum dot formed in a Si/SiGe 
heterostructure. We present measurements of transport 
through the device, demonstrating controllable tunnel 
coupling between the quantum dot and the localized 
state. A set of stability diagram measurements enable 
a determination of the relative magnitude of the capac¬ 
itance between the surface gates and both the quantum 
dot and the localized state. We report the expected elec¬ 
tron density profiles in the quantum dot and the neigh¬ 
boring reservoirs. Combining the experimental results 
with 3D capacitive modeling based on the electron den¬ 
sity profiles, we determine the most likely location of the 
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localized state in the device. The results reported here 
demonstrate that it is possible to control the tunnel rate 
between localized states and quantum dots, even though 
there is a dramatic difference in the characteristic length 
scale between these two electronic systems. 

A gate-defined quantum dot was fabricated in a 
Si/Sio. 68 Geo .32 heterostructure grown by chemical vapor 
deposition on a relaxed buffer layer in which the surface 
roughness had been removed by chemical-mechanical pol¬ 
ishing. A schematic diagram of a cross-section of the de¬ 
vice is shown in Fig. 1(a), and a top-view, false-color 
scanning electron micrograph of the sample is shown 
in Fig. 1(b). Measurements were performed in a di¬ 
lution refrigerator with a mixing chamber temperature 
T M c < 30 mK. 

Figure 1(c) is a Coulomb diamond plot of the cur¬ 
rent through the device: it shows the derivative of the 
source-drain current between ohmic contacts S and D, 
as labelled in Fig. 1(b), as a function of the voltage Vq 2 
on gate G2. The primary features in the plot are the 
diamond-shaped regions of very low differential conduc¬ 
tance characteristic of Coulomb blockade with an av¬ 
erage charging energy E c = e 2 /CDot = 760 peV. By 
comparing this charging energy with the charging ener¬ 
gies of few-electron Si/SiGe quantum dots 24 , it is clear 
that this quantum dot is in the many electron regime. 
We can also obtain, from the excited states visible in 
Fig. 1(c), an estimate of the single-particle energy of 
about 380 peV in this quantum dot. The data also en¬ 
able the extraction of the proportionality constant (the 
lever arm) ac 2 = 148 peV/mV between the voltage on 
gate G2 and the energy of the quantum dot. 

In addition to these expected features, the red arrows 
near the center of Fig. 1(c) highlight a sharp, isolated 
charging event. Because this feature is so isolated—over a 
range in gate voltage corresponding to the addition of 13 
electrons to the main dot only this one additional feature 
is observed—it corresponds to a very small, large charg- 
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FIG. 1. Device design and characterization data of a quantum dot, revealing evidence of a nearby localized state, (a) Schematic 
diagram showing a side view of the device. Substrate is a Si/Sio. 6 sGeo .32 heterostructure with a 10 nm silicon well (light grey) 
and 32 nm SiGe offset (dark grey). Both the upper (purple and red) and lower (green) layers of gates are 2 nm titanium and 
20 nm gold deposited by electron beam evaporation. The lower gates were deposited on 10 nm of atomic layer depostion (ALD) 
grown aluminum oxide (light orange) while the upper gates were on 90 nm of oxide. Ohmic contacts S and D (denoted with M 
symbols) are 5 nm titanium and 40 nm gold on a region of the heterostructure degenerately doped with phosphorus through 
the quantum well (black dashed boxes). Approximate location of the quantum dot and impurity are shown schematically in 
the center of the figure (yellow dashed oval), (b) Scanning electron micrograph of a completed device identical to the measured 
device. Upper gate A (shown in purple) and paddle gate P (shown in red) were positively biased to accumulate a 2DEG in 
the reservoir and to control the energy of the dot, respectively. On the lower level, gates Gl, G2, and QPC (shown in yellow) 
were negatively biased to provide the dot confinement potential, whereas gates B1 and B2 (shown in green) controlled the 
tunnel barriers to the source (S) and drain (D) ohmic contacts (denoted with M symbols), (c) The derivative dIsr>/dVsr> of the 
transport current with respect to the gate voltage Ug 2 - Multiple Coulomb diamonds are observed. Near the center of the plot, 
a sharp resonance (indicated by arrows) is observed, suggesting a localized state, (d) The current /sd at fixed Vsd — 100 pV 
as a function of Up and Ug 2 , showing many charge transitions of the dot. A jump in gate voltage of the location of the Coulomb 
blockade transitions is observed (indicated by arrows), corresponding to the localized state in (c). 


ing energy object. In addition, the capacitance between 
various gates and the object corresponding to this fea¬ 
ture are different from those corresponding to the quan¬ 
tum dot. This difference is made immediately clear in 
Fig. 1(d), which reports the current through the quan¬ 
tum device as a function of the voltages on gate P and 
G 2 . The phenomenology of this plot is very similar to 
those observed in metal-oxide-semiconductor devices in 
which donors have been implanted : 25 near the center of 
the scan, a series of shifts in the charge transitions of 
the dot can be observed; these shifts correspond to the 
feature marked by the red arrows in Fig. 1(c). The line 
through the gate voltage space spanned by Vp and Vq 2 
connecting these shifts has a different slope than that of 
the Coulomb blockade peaks corresponding to the dot, 
confirming the presence of a nearby localized state that 
is not at the same physical position as the dot. 

Figure 2(a) reports a high-resolution measurement of 
the source-drain current /sd as a function of Up and 
Vbi- The dashed black lines indicate the voltages at 
which the localized state charges, corresponding to the 
observed shift and gap in the Coulomb blockage peaks 
corresponding to the quantum dot. No current is ob¬ 
served along these black dashed lines, indicating that— 
unlike the case for the quantum dot—the localized state 
is not tunnel coupled to both the source and the drain. It 
is possible, however, that the localized state is connected 
to either the source or the drain, and this hypothesis is 
supported by the faintly visible line of current (dark blue 


in the color scale) that sits at the position of the polar¬ 
ization line in this stability diagram. 

Figs. 2(b) & (c) show very high resolution measure¬ 
ments of a pair of triple points in this two-site system, 
for two different values of the voltage on gate B2. In both 
plots, there is no current along the black dashed lines cor¬ 
responding to the charge transition of the localized site, 
confirming that the localized site is not tunnel-coupled 
to both the source and the drain. However, in Fig. 2(b), 
where Up 2 = —402 mV, current is observed along the 
polarization line; in contrast, no such current is observed 
in Fig. 2(c), where Ub 2 = —404 mV. This current is 
studied in more detail in Fig. 2(d), where we report line 
cuts across the polarization line, as indicated by the gray 
dashed line in Fig. 2 (b). As a function of gate voltage 
Up 2 5 Fig. 2(d) shows a dramatic evolution of the current 
along this path. For Vp 2 equal to either -398 or -401 mV, 
no peak in current occurs at the polarization line. In con¬ 
trast, for intermediate values of Vp 2 = -399 and -400 mV, 
a prominent peak in current is observed at the position 
of the polarization line. The narrow voltage range of the 
barrier gate over which the tunneling condition is satis¬ 
fied suggests that transport through the impurity is very 
sensitive to the relative strengths of the tunnel couplings 
of the dot and the impurity to each other as well as to 
the leads. 

The difference between the current peak shown in 
Fig. 2(d) and the conventional Coulomb peaks corre¬ 
sponding to the quantum dot is also highlighted by 
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FIG. 2. Data demonstrating control of the tunnel coupling 
between the localized state and the quantum dot. (a) A high- 
resolution data set showing jumps in the Coulomb blockade 
transitions as a function of Vp and Vbi. A double dot-like sta¬ 
bility diagram is revealed, including a slightly visible current 
on the polarization line. The single electron transitions of the 
localized state are not visible in the current /sd and instead 
are shown schematically as black dashed lines, (b) &; (c) 
Very high resolution scans of one of the state crossings shown 
in (a), for two different values of the gate voltage Vb 2 • In (b), 
Vb 2 = —402 mV and current is present at the polarization 
line. In (c), Vb 2 — —404 mV, and current is absent at the 
polarization line, (d) Line scans along a path shown by the 
light gray line in panel (b). The data for the Vb 2 — —398 mV 
trace are not shifted, and subsequent traces are offset ver¬ 
tically by 100 f A and laterally by 1 mV each for visibility. 
These plots, which were acquired in close succession in time, 
show that changing Vb 2 changes the tunnel couplings to the 
localized state, turning on and off current at the polarization 
line. (The overall conditions in this plot are slightly different 
than those in panels b and c.) (e) Black inverted triangles 
show the height of the main Coulomb blockade peak near the 
anticrossing with the localized state, and red triangles show 
the height of the current peak on the polarization line, both of 
which are strongly temperature dependent. Inset: Coulomb 
blockade peak height Ip far from the anticrossing with the lo¬ 
calized state. The temperature dependence in this regime is 
conventional for a large dot at reasonably low temperature. 26 


the temperature dependence of each peak. The inset 
to Fig. 2(e) shows the temperature dependence of the 
Coulomb blockade peak for the single dot, which is tun¬ 
nel coupled to both the source and the drain, for gate 
voltages such that the localized state is not involved in 
the transport. The current is nearly constant as a func¬ 
tion of temperature, rising slightly as the temperature 
drops, consistent with a reasonably large quantum dot 
at temperatures T for which kT E c (Ref. 26). In 


the main panel of Fig. 2(e), the black inverted triangles 
show the temperature dependence of the main Coulomb 
blockade peak very close to the anticrossing with the lo¬ 
calized state, and the red triangles show the temperature 
dependence of the current peak on the polarization line. 
In contrast with the behavior in the inset, both of these 
peaks increase strongly with increasing temperature, con¬ 
sistent with activated behavior. This behavior is consis¬ 
tent with a localized state tunnel coupled to the dot and 
one (but not both) of the reservoirs. Considering first the 
main Coulomb peak: on this peak, transport through the 
dot is allowed (by definition), but accessing the localized 
state requires thermal activation, which when present at 
elevated temperatures opens a parallel path to exit the 
dot, increasing the total current. Considering the polar¬ 
ization line peak (which is more than 5 times weaker than 
the main Coulomb peak): this current is suppressed at 
low temperature. Although charge can shuttle between 
the dot and the localized state at no energy cost, it can¬ 
not tunnel to or from the leads—raising the temperature 
activates this process. 


To better understand these results, we combined sys¬ 
tematic experimental stability diagram measurements 
with electrostatic device modeling to determine the loca¬ 
tion of the localized state 27 . Experimentally, we defined a 
central operating point in gate voltage space near the lo¬ 
calized state-to-quantum dot charge transition. Around 
this point, as shown in Fig. 3(a), we performed five two- 
dimensional stability diagram measurements, sweeping 
gate voltage Vp, which we use as our reference, as we 
stepped five other voltages: Vq i, Vbi, Vg2, Vqpc , and 
Vb 2 - For each scan, all other voltages were held fixed to 
their values at the central operating point. 

To interpret these data, we constructed an electro¬ 
static device model in COMSOL Multiphysics 28 , using 
the device geometry from the experiment, as illustrated 
in Fig. 3(b). The model was solved in the Thomas-Fermi 
Approximation 29,30 , with the self-consistent charge ac¬ 
cumulation determined using a 2D density of states con¬ 
fined to a sheet at the Si-SiGe interface. We assumed a 
2x2 carrier degeneracy due to (2) spin and (2) valley 
degrees of freedom. The computational domain was 3x5 
fi m laterally, and included a 100 nm air cap above the 
oxide layer and 500 nm SiGe substrate below the silicon 
well. We used zero-field boundary conditions on the sides 
of the domain and the top of the air cap, and set the con¬ 
duction band edge to the Fermi level at the bottom of the 
domain. In addition, we used voltages: V G 2 = —0.335 V, 
Vb2 = -0.400 V, Vqpc = -0.100 V, V G1 = -0.104 V, 
Vbi = +0.270 V, V A = +2.75 V, and V P = +2.75 V. 

We approximate the dot and 2DEG regions predicted 
by COMSOL as 5 nm thick metallic sheets at the 5 x 10 11 
cm -2 density contour (Fig. 3(b)), and we treat the lo¬ 
calized state as a 1 nm radius metallic sphere. For a 
given placement of the localized state, we construct a ca¬ 
pacitance model that predicts each of the experimental 
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FIG. 3. Locating the localized state by combining experimen¬ 
tal data and an electrostatic device model, (a) A series of 
two-dimensional gate scans at fixed Vs d measuring transport 
through the quantum dot-impurity stability diagram. Each 
plot reports the result of a different gate voltage being swept 
on the horizontal axis, versus the gate voltage Vp on the verti¬ 
cal axis. The black dashed lines highlight the Coulomb block¬ 
ade transitions of the dot while the white dashed lines show 
the expected position of the unseen localized state charging 
event, (b) Results of 3D electrostatic modeling of the exper¬ 
iment to determine the location of the localized state. The 
top layer shows the gate geometry of the device local to the 
dot. The middle layer shows the electron densities of the dot 
and leads as calculated by COMSOL, with the contour corre¬ 
sponding to 5 x 10 11 cm -2 electron density. The bottom layer 
of the device shows the most likely location of the impurity, as 
determined by the discrepancy metric (Eq. (1)), directly un¬ 
der the tip of gate G1 near the lower Si/SiGe interface about 
10 nm below the dot and 2DEG. 


charge-stability diagrams 31 . By rastering the localized 
state position across the device and computing as a fit 
metric a weighted sum-of-squared-differences between ex¬ 
perimental and predicted values, we estimate the location 
of the localized state that is most consistent with the ob¬ 
served data shown in Fig. 3(a). In particular, we sum the 


squared differences of two types of observed quantities: 

1. The slope of the line connecting all of the “offsets” 
in the Coulomb blockade lines (the white dashed 
lines in Fig. 3(b)). 

2. The magnitude of the jump along the y-axis of a 
Coulomb blockade line due to the localized state. 

Quantities of type 1 are unitless slopes whereas quantities 
of type 2 have units of energy. To combine these into a 
single discrepancy metric, we found empirically that we 
needed to scale the type 2 quantities reported in meV 
by 1 x 10 8 to balance them with the quantities of type 
1. Thus, the overall discrepancy metric that is plotted in 
Fig. 3(b) is equal to: 

5 5 

DIS = Y J { p Oi-EO i ) 2 +lxl^xY J {PJi-E J i ) 2 , (1) 

i=1 i=1 

where i indexes each of the five experimental slices shown 
in Fig. 3(a), PO and EO are the predicted and experi¬ 
mental type 1 offset quantities, and PJ and EJ are the 
predicted and experimental type 2 jump quantities re¬ 
spectively. The lower layer of Fig. 3(b) shows a cut of 
DIS along a plane 12 nm beneath the top of the strained 
Si well, identifying a region under the tip of gate G1 as 
the most likely region in the x-y plane to find the local¬ 
ized state. We find that the DIS value is not very sensi¬ 
tive to the depth (z-coordinate) between 10 and 20 nm, 
and thus we show the 12 nm data since this appears most 
consistent with the observed tunneling through the state. 

Considering measurements and simulations together, 
we propose a tunnel rate dependent model of our hybrid 
quantum dot-impurity system. Under typical device op¬ 
eration, the impurity is tunnel coupled to one of the leads 
and capacitively coupled to the dot. Coupled this way, 
a charging event of the localized state varies the electric 
field local to the dot, changing the Coulomb blockade 
condition and resulting in the familiar jump in the dot 
charge transition from Fig. 2(c). Changing the voltage of 
the tunnel barrier Vb 2 changes the dot-drain, impurity- 
drain, and dot-impurity tunnel rates. Under certain tun¬ 
ings, like those shown in Fig. 2(b), the three tunnel rates 
allow for current through the normally blockaded region 
as well as enhancement of current corresponding to the 
dot charge transitions. 

In conclusion, we have shown measurements and mod¬ 
eling of a tunnel coupled quantum dot-impurity system 
in a Si/SiGe heterostructure. We demonstrated tunable 
tunnel coupling between the impurity and the dot that 
is controlled by varying a nearby gate voltage, and we 
reported the temperature dependence of the coupled sys¬ 
tem. We also have found the most likely position of the 
localized stated through capacitive modeling, and the ca¬ 
pacitances extracted from this model are in good agree¬ 
ment with the experimental results. 
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